Boundary problems for Dirac electrons and edge-assisted Raman scattering in 

graphene 



D. M. BaskcEl 

Laboratoire de Physique et Modelisation des Mileux Condenses, 
Universite Joseph Fourier and CNRS, 25 Rue des Martyrs, BP 166, 38042 Grenoble, France 

The paper reports a theoretical study of scattering of electrons by edges in graphene and its effect 
on Raman scattering. First, we discuss effective models for translationally invariant and rough 
edges. Second, we employ these models in the calculation of the edge-activated Raman D peak 
intensity and its dependence on the polarization of the incident/scattered light, as well as on the 
position of the excitation spot. Manifestations of the quasiclassical character of electron motion in 
Raman scattering are discussed. 



I. INTRODUCTION 

Graphene, a monolayer of carbon atoms, first ob- 
tained in 2004fi received increasing interest due to its 
exceptional electronic properties and original transport 
physics,^ Gradual miniaturization of graphene devices in- 
creases the importance of edge effects with respect to the 
bulk two-dimensional physics. Starting from a graphene 
sheet, nanoribbons and quantum dots can be produced 
by lithography and etching.— i^i^i^i^ii^ Edges also play a 
fundamental role in the quantum Hall effect. 

Graphene edges can be studied by several experimental 
techniques. Scanning tunneling microscopy (STM) and 
transmission elecron microscopy (TEM) can resolve the 
structure of the edge on the atomic scalci^ i^^'^^i-'^^'^'^i^^ 
Raman scattering has also proven to be a powerful tech- 
nique to probe graphene edges i^^i^^'^^i"'^^'^^ The so-called 
D peak at 1350 cm~^ is forbidden by momentum con- 
servation in a perfect infinite graphene crystal, and can 
only be activated by impurities or edges. Invoking the 
double-resonance mechanism for the D peak activation;^ 
Cangado et al. have shown that a perfect zigzag edge 
does not give rise to the D peakji^ It should be empha- 
sized that this property is determined by the effect of the 
edge on the electronic states. 

A great deal of theoretical studies of electronic proper- 
ties near the edge has focused on the case of ideal zigzag 
or armchair edges, most commonly adopting the tight- 
binding description. One of the spectacular results ob- 
tained by this approach was the existence of electronic 
states confined to the zigzag edgej^i^i^^i^ii^^ which 
was later confirmed exp er iment allyi^ ' i " 1 1 1 1 1 The question 
about general boundary condition for Dirac electron wave 
function at a translationally invariant graphene edge has 
been addressed^ and a detailed analysis of boundary 
conditions which can arise in the tight-binding model has 
been performedi^ In spite of the fact that all graphene 
samples produced so far have rough edges, the number 
of theoretical works dedicated to rough edges is limited. 
Most of them model edge roughness in the tight-binding 
model by randomly removing lattice sitesi^^i^i^Si^ The 
opposit limit of smooth and weak roughness has been 
consideredi^ Edge states on zigzag segments of finite 
length have also been studied recentlyi^^ 



The present work has several purposes. One is to de- 
velop analytically treatable models which would describe 
electron scattering on various types of edges in terms of 
as few parameters as possible. The second one is to calcu- 
late the polarization dependence of the D peak intensity 
for different models of the edge, and thus see what infor- 
mation can be extracted from this dependence. The third 
one is to identify the characteristic length scale which 
confines the Raman process to the vicinity of the edge, 
i. e. the spatial extent of the Raman process. It will be 
shown that the last two issues are intimately related to 
the quasiclassical character of the electron motion during 
the Raman scattering process. 

The paper is organized as follows. In Sec. |TT] we dis- 
cuss the problem in qualitative terms and summarize 
the main results of the work. In Sec. IIIII we summarize 
the Dirac description of single-electron states in an infi- 
nite graphene crystal and formulate the Huygens-Fresnel 
principle for Dirac electrons. In Sec. II VI we discuss mod- 
els for the electron scattering from a graphene edge, con- 
sidering translationally invariant as well as rough edges. 
Sec.lVlintroduccs the model for electron-phonon coupling 
and describes the general scheme of the calculation of the 
D peak intensity using the standard perturbation theory 
in the coordinate representation. Finally, Sees. IVII IVIII 
and IVIIII are dedicated to the calculation of the D peak 
intensity for an ideal armchair edge, an atomically rough 
edge, and an edge consisting of a random collection of 
long zigzag and armchair segments, respectively. 



II. QUALITATIVE DISCUSSION AND 
SUMMARY OF THE MAIN RESULTS 

A. Electron scattering by the edge 

First, we discuss translationally invariant edges. For 
example (see Fig. [T]), a zigzag edge has a spatial pe- 
riod de = aVS (a K. 1.42 A is the C-C bond length), 
an armchair edge has de ~ 3a, and a more compli- 
cated edge, shown in Fig. [Ijc), has de = V^a « 4.6a 
(the spatial period is measured along the average di- 
rection of the edge). It is important to compare de to 
the electronic wavelength (we prefer to divide the lat- 
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FIG. 1: Examples of ordered edges: (a) zigzag, (b) armchair, 
(c) a more complicated but still translationally invariant edge. 



valleys are mixed randomly, as there is no small or large 
parameter which would suppress or favor any particular 
channel (the smallness a/A^ ^ 1 suppresses the direct 
intravalley scattering, but multiple intervalley scattering 
efficiently mixes states within the same valley as well). 
The electron is thus scattered randomly both in all di- 
rections (as a consequence of the randomization of the 
momentum direction within the same valley), and be- 
tween the two valleys. For this case in Sec. IIVBI we 
propose a phenomenological model which describes such 
random scattering of electrons by the edge, respecting 
only the particle conservation and the time-reversal sym- 
metry. Essentially, each point of the edge is represented 
by an independent point scatterer which randomly ro- 
tates the valley state of the electron. This model is used 
for quantitative calculations in the subsequent sections. 

Edges, rough on length scales much larger than the 
lattice constant, are likely to consist of distinct segments 
of zigzag and armchair edges, as shown by STMiSiiiii^ 
and TEMJ^ii^ Then the overall probability of scatter- 
ing within the same valley or into the other valley is 
simply determined by the fraction of the corresponding 
segments. The problem of angular distribution of the 
scattered electrons is analogous to the well-studied prob- 
lem of light scattering by rough surfacesj ^"^'^^ '^'^'^'^'^ 
The main qualitative features of the scattering, namely, 
the sharp coherent peak in the specular direction, the 
smooth diffuse background, and the enhanced backscat- 
tering peak, should be analogous for the electrons in 
graphene as well. The so-called surface polaritons, shown 
to play an important role in the light scattering, are anal- 
ogous to the electronic edge states in graphene. Still, full 
adaptation of this theory for the case of Dirac electrons 
in graphene represents a separate problem and is beyond 
the scope of the present work. Here we only consider the 
case when regular edge segments are sufHciently long, i. e. 
their typical length Xe- Then the diffraction correc- 
tions, small in the parameter X^/de ^ 1, can be simply 
found as it is done in the classical optics;^ using the 
Huygens-Fresnel principle for Dirac electrons, Eq. (j3.7p . 



ter by 27r), Ag = f/|e|, where e is the electron energy, 
and V ~ l.l ■ 10* cm/s « 7.3 eV • A is the electron ve- 
locity (the slope of the Dirac cones). For comparison, 
at e = 1 eV Ag « 7.3 A « 5a. As long as de < ttAc, 
the component of the electronic momentum along the 
edge, p|| , is conserved (we measure the electron momen- 
tum from the Dirac point). For longer periods the edge 
acts analogously to a reflective diffraction grating in op- 
tics; this case is not considered here. In the limit d^ <^ Ae 
the reflection of electrons from any periodic edge can be 
described by an effective energy-independent boundary 
condition for the electronic wave functioni^i21 

Next, we study rough edges. An extreme case is when 
the edge is rough at the atomic scale, like that in the 
tight-binding model with randomly removed sites. Then 
it is reasonable to assume that in the vicinity of the edge 
all plane-wave components with a given energy in both 



B. Quasiclassical picture of Raman scattering 

Since graphene is a non-polar crystal, Raman scatter- 
ing involves electronic excitations as intermediate states: 
the electromagnetic field of the incident laser beam inter- 
acts primarily with the electronic subsystem, and emis- 
sion of phonons occurs due to electron-phonon interac- 
tion. The matrix element of the one-phonon Ramanpro- 
cess can be schematically represented as 

J. {i\He-em\a){a\He-ph\b){b\He-er,i\f) 

Here \i) is the initial state of the process (the incident 
photon with a given frequency and polarization, and no 
excitations in the crystal), |/) is the final state (the scat- 
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tered photon and a phonon left in the crystal), while 
I a) and \b) are the intermediate states where no photons 
are present, but an electron- hole pair is created in the 
crystal and zero phonons or one phonon have been emit- 
ted, respectively. Note that these intermediate states 
correspond to electronic eigenstates in the presence of 
the edge, i. .e, scattered states rather than plane waves. 
Ei = Ef, Ea, and Et are the energies of the correspond- 
ing states, and 27 is inverse inelastic scattering time (the 
overall rate of phonon emission and electron-electron col- 
lisions). He-em and H^-ph stand for the terms in the sys- 
tem hamiltonian describing interaction of electrons with 
the electromagnetic field and with phonons, respectively. 

As discussed in Refs. li^lisl . for one-phonon scattering 
processes it is impossible to satisfy the energy conser- 
vation in all elementary processes. This means that the 
electron-hole pair represents a virtual intermediate state, 
and no real populations are produced. Formally, at least 
one of the denominators in Eq. ()2.1|1 must be at least 
of the order of the phonon frequency Wph 3> 7. In fact, 
the main contribution to the matrix element comes from 
such states that the electron and the hole have the en- 
ergy e close (within ~ Wph) to half of the energy ujin of 
the incident photon: |e — ti;m/2| ~ Wph. These two energy 
scales are well separated: Wph « 0.17 eV, while typically 
ijJin/2 « 1 eV. According to the uncertainty principle, the 
energy uncertainty, Wph, determines the typical lifetime 
of the virtual state (electron- hole pair), ^ 1/wph- This 
time scale determines the duration of the whole process. 

As we are dealing with a translationally non-invariant 
system, it is useful to analyze the Raman process in the 
coordinate representation. The time scale l/wph, intro- 
duced above, translates into the length scale £ph = v/uiph- 
Thus, this length scale, ^ph ~ 4 nm, determines the spa- 
tial extent of the process (we will return to this point be- 
low). Its largness compared to the electron wavelength, 
^ph/^e = ij-'in/(2wph) ^ 1, cnsurcs that the electronic 
wave functions determining the matrix elements for each 
elementary process, admit a quasiclassical representa- 
tion. The quasiclassical approximation for the electronic 
wave functions is fully analogous to the geometrical op- 
tics approximation for electromagnetic waves, electronic 
trajectories corresponding to light rays. Corrections to 
this approximation are known as diffraction and are small 
by the parameter ujph/ojin <C 1. It should be emphasized 
that the quasiclassical picture is neither an assumption, 
nor a hypothesis, but it arises automatically in the di- 
rect calculation of the Raman matrix element which is 
performed in the main part of the paper. 

In the quasiclassical picture, the photoexcited electron 
and hole can be viewed as wave packets of the size ~ , 
initially created at an arbitrary point of the sample. More 
precisely, instead of a point one can consider a region of a 
size 61, such that ^ 61 <ti £ph- Then momentum con- 
servation holds up to (5p ~ 1/61 ^ e/v by virtue of the 
uncertainty principle, so that electron and hole momenta 
whose magnitude is e/v (counted from the Dirac point), 
have approximately opposite directions, as the photon 




FIG. 2: (Color on-line) Real-space representation of the scat- 
tering process responsible for the D-peak near graphene edges. 
The lightning represents the incoming photon which gener- 
ates the electron-hole pair. The solid black arrows repre- 
sent the quasi-classical trajectories of the electron and the 
hole. The dashed arrow represents the emitted phonon. The 
flash represents the radiative recombination of the electron- 
hole pair producing the scattered photon, (a) Backscattering 
off a translationally invariant edge is possible only at normal 
incidence (up to the quantum uncertainty), (b) For oblique 
incidence on a translationally invariant edge the reflection is 
specular, so the electron and the hole will not be able to meet 
at the same point, (c) For a rough edge backscattering is 
possible even at oblique incidence. 



momentum is very small. The same argument holds for 
the phonon emission and for the radiative recombination 
process: in order to emit a photon, the electron and the 
hole must meet in the same region of space of the size 61 
with almost opposite momenta (up to 1/61). Momentum 
conservation at the reflection from the edge depends on 
the quality of the edge, as discussed in Sec. Ill Al Regard- 
less of the properties of the edge, an elementary geomet- 
ric consideration, illustrated by Fig. [2l shows that for the 
electron and the hole to be able to meet, the scattering 
on both the phonon and on the edge must be backward. 

In the quasiclassical picture, the electron and the hole 
have to travel the same distance between creation and 
annihilation, as their velocities are equal. Then the pro- 
cess in Fig. [3] (a) has more phase space satisfying this 
restriction, and this process gives the main contribution 
to the Raman matrix element. This will be also shown 
by an explicit estimate in Sec. IVIBI Note that the three 
processes shown in Fig. [31 can be considered in the mo- 
mentum space, as shown in Fig. [4] According to the 
abovesaid, the processes (b) and (c), often shown in the 
literature as an illustration of the double resonance,— are 
in fact weaker than the process (a) by a factor ~ Wph/i^m- 



4 




FIG. 3: (Color on-line.) Real space representation of different 
contributions to the matrix element of the scattering process 
responsible for the D-peak at an ideal armchair edge, placed 
at a; = 0. The solid black arrows represent the quasi-classical 
trajectories of the electron and the hole corresponding to the 
three Green's functions in Eqs. (|5.7c|l . (|5.7d[l . Trajectories 
(a), (b), (c) correspond to decomposition of each of the three 
Green's functions in Eq. (|5.7c|) . t is the overall spatial extent 
of the process, and Ae = v/e is the electron wavelength divided 
by 2-K. 




FIG. 4: (Color on-line.) Momentum space representation of 
different contributions to the matrix element of the scattering 
process responsible for the D-peak at an armchair graphene 
edge, corresponding to the real space picture shown in Fig. O 
Solid lines represent the Dirac cones around K and K' points 
of the first Brillouin zone. Vertical solid arrows represent 
interband electronic transitions accompanied by photon ab- 
sorption or emission (photon wave vector is neglected) , dashed 
arrows represent phonon emission, the horisontal dotted ar- 
row represents the scattering from the edge. 



C. Polarization dependence 

The backscattering condition has immediate conse- 
quences for the polarization dependence of the Raman 
scattering intensity. Indeed, the matrix element of cre- 
ation/annihilation of an electron and a hole pair with 
momenta p, — p (counted from the Dirac point) by a pho- 
ton with the polarization e, is proportional to [e x p]^, 
reaching its maximum when e _L p. Since a perfect edge 
conserves the component of momentum along the edge, 
the backscattering is possible only at normal incidence, 
as seen from Fig. [2]).— This gives the polarization depen- 
dence of the D peak intensity as Id oc sin^ ipin siv? ^pout , 
where and ipout sue the angles between the polar- 
izations of the incident and scattered photons and the 
normal to the edge. If one does not use the analyzer 
to fix the polarization of the scattered photons, the de- 
pendence is Id oc sin^ ipin- In experiments, however, the 
intensity never goes exactly to zero for the polarizations 
perpendicular to the edge, but remains a finite fraction e 
of the intensity in the maximum.— li^ii^ What determines 
this residual intensity in the minimum? 

For an ideal edge the finite value of the intensity in 
the minimum is entirely due to the quantum uncertainty. 
Namely, the momenta of the electron and the hole upon 
their creation are not exactly opposite, but up to an 
uncertainty ~ 1/^'; the annihilation occurs not exactly 
at the same spatial point, but within the spatial uncer- 
tainty 61. If the spatial extent of the process is £, the 
uncertaintly is estimated as 61 \/A7?, and the ratio 
e of the intensities in the minimum and in the maxi- 
mum (i. e., for polarizations perpendicular and parallel 
to the edge, respectively) should be small as a power of 
the small parameter quasiclassical parameter X^/£. The 
calculation is performed in Sec. IVIBl the result is given 
by Eqs. (|6.5b[) and (|6.5cp for the detection without and 
with the analyzer, respectively. Up to logaritmic factors, 
the ratio e ~ i^ph/^in ~ (^e/^ph)^- This corresponds to 
£ph ~ £, in accordance with the energy-time uncertainty 
principle, as discussed in the previous subsection. 

For rough edges the intensity in the minimum is deter- 
mined by the ability of the edge to backscatter electrons 
at oblique incidence, as shown in Fig. [UJc). If the edge 
is rough at the atomic scale, oblique backscattering is 
nearly as efficient as normal backscattering. Still, such 
oblique trajectories are longer than those corresponding 
to the normal incidence, so they are expected to have a 
smaller weight since the virtual electron-hole pair lives 
only for a restricted time. So one still can expect a mini- 
mum of intensity for perpendicular polarization, but it is 
of a purely geometrical origin, so one does not expect a 
parametric smallness of the ratio s. The calculation using 
the model Sec. lIVBl for an atomically rough, is performed 
in Sec. IVII K[ and the result is given by Eqs. (j7.3bp and 
(|7.3cP for the detection without and with the analyzer, 
respectively. In the former case the ratio s ~ 1/3, and 
the minimum is indeed for the polarization of the inci- 
dent light perpendicular to the edge. With the analyzer. 
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the absolute minimum is e = 1/10, reached when the 
polarizer and the analyzer are oriented at the angle 7r/3 
with respect to the edge and to each other. When the 
polarizer and the analyzer arc both rotated parallel to 
each other, the minimum is e = 1/5. 

For an edge consisiting of segments longer than elec- 
tronic wavelength, dg S> Ac, we first analyze the contri- 
bution of a single armchair segment. It is calculated in 
Sec. IVIIIl and given by Eqs. ((53a|) and ((83b| for the 
detection without and with the analyzer, respectively. 
The minimum is reached for the polarization, perpendic- 
ular to the armchair direction (which does not have to 
coincide with the average direction of the edge!), and is 
determined by the quantum diffraction of the electron on 
the segment, e ~ Xc/de (provided that X^/de > cjpj^/wf^, 
the ratio for the infinite armchair edge, which is the case 
when de ^ 50 nm for ujm — 2eV). To obtain contribution 
of the whole edge, it is sufficient to multiply these expres- 
sion by the total number of such segments and replace d^ 
by its average, provided that all armchair segments have 
the same orientation. 

It is crucial, however, that up to three different orien- 
tations of armchair segments are possible, at the angle 
7r/3 to each other. When contributions corresponding 
to several different orientations are added together, the 
polarization dependence may change quite dramatically, 
as was briefly discussed by the author and coworkers in 
Ref. [l^ and is considered in more detail in Sec. IVIIIl Note 
that if the average direction of the edge is armchair or 
zigzag, the possible orientations of armchair segments are 
symmetric with respect to the average direction: three 
orientetaions at angles and ±7r/3 for an armchair edge, 
and two at angles ±7r/6 for a zigzag edge. Most likely, the 
number of segments corresponding to "-f " and "— " signs 
will be equal on the average. Then, by symmetry, the 
maximum of intensity in both cases will be reached for 
the polarization along the average direction of the edge, 
in agreement with recent experimental observations . 

When the average direction is armchair, the ratio be- 
tween the minimum and the maximum of intensity is de- 
termined by the relative fraction of segments oriented at 
±7r/3 with respect to those oriented along the average 
direction. When the average direction is zigzag, the po- 
larization dependence is fully determined by the symme- 
try (on the average) between the two armchair directions, 
and is given by Eq. (|8.10p . The ratio between the mini- 
mum and the maximum is e = 1 /3 for detection without 
analyzer, and e — 1/9 for detection with an analyzer 
parallel to the polarizer of the incident light, again, in 
agreement with Refs. 18.19. Thus, quantum diffraction 
effects appear to be masked by the purely geometrical 
effects. 

Remarkably, the quantum diffraction limit is still ac- 
cessible if only two orientations of armchair segments are 
present (which is the case when the average direction is 
zigzag or close to it). It is sufficient to put the polarizer 
perpendicular to one of the armchair directions, and the 
analyzer perpendicular to the other one, thereby killing 



the leading specular contribution for both segments. In 
this polarization configuration the absolute minimum of 
the intensity is reached, and it is indeed determined by 
the quantum diffraction, as given by Eq. (jS.llap . 



D. Excitation position dependence 

In Ref.lH, confocal Raman spectroscopy was suggested 
and used as a way to probe the length scale £ which re- 
stricts the D peak to be in the vicinity of the edge (the 
spatial extent of the Raman process). The idea is to focus 
the incident light beam as tightly as possible, so that its 
electric field £i„ (r) has the shape 5i„ (r) cx e"''""''''! /^^^ \ 
where the width L can be measured independently, and 
the spot center position Tq can be varied experimentally. 
Then, given that the intensity of the D peak is propor- 
tional to 

/i3cx J IC{r,r')£Ur)£*n{r')d^rd\', (2.2) 

where /C(r, r') is a certain kernel, decaying away from the 
edge with a characteristic length scale £, a measurement 
of the dependence /^(ro) would give information on the 
kernel. The measurement would be especialy simple in 
the case L <^ i. In reality, however, the relation is the 
opposite: in Ref. [l^ L = 186.5 nm, and ^ is a few tens of 
nanometers at most 4^ 

In this situation the dependence of the Raman inten- 
sity Id{yo) is very close to the excitation intensity profile 
|fi„(r)p, and the nonlocality of the kernel /C(r,r') mani- 
fests itself only in a slight change of shape of/D(ro) with 
respect to |fi„(r)p. In the first approximation it can 
be viewed just as small shift and broadening. When the 
signal-to-noise ratio is not sufficiently high to perform the 
full functional deconvolution, one has to assume a specific 
functional form for the kernel and do a few-parameter fit. 
It is clear that different functional forms will give values 
of £, differing by a factor of the order of 1. In Ref. [l^ 
the form /C(r, r') — 9{x) 9{x') was assumed, 

where x is the distance from the edge, 0{x) is the step 
function, and £^ — v/{2j) is the electron inelastic scat- 
tering length [27 is the electron inelastic scattering rate, 
see Eq. (j2.ip ]. This assumption seems to contradict the 
fact that the lifetime of the virtual electron-hole pair is 
^ 1/ ojoh, d iscussed in Sec. Ill B[ as it was pointed out in 
Refs. ^MM 

The explicit form of the kernel /C(r, r') for an ideal 
armchair edge is calculated in Sec. IVI CI it is given by 
Eq. (j6.7cp . and it turns out to be more complicated than 
a simple exponential. In fact, it depends on both length 
scales, £p\i and £^. The length is shorter, but the 
spatial cutoff it provides is only power-law. The longer 
length £ry is responsible for the strong exponential cut- 
off. Which of the two lengths plays the dominant role 
in the Raman process, turns out to depend on the spe- 
cific observable to be measured. The total integrated 
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intensity for 5m (r) = const is proportional to the in- 
tegral J IC{r,r') cPr (Pr' , which is determined mainly by 
^ph, while £j enters only in a logarithmic factor. The 
same can be said about the polarization dependence and 
diffraction corrections, discussed in the previous subsec- 
tion. However, the change in the shape of /^(ro), com- 
pared to the excitation intensity profile |£i„(r)p, is deter- 
mined by the second and higher moments of the kernel, 
/ a;"(x')"' JC{r,r')(fr(fr', with n + n' > 2. These mo- 
ments turn out to be determined by the longer scale 
Thus, the interpretation by the authors of Ref. [TH of their 
experiment is qualitatively correct. 

Analysis of the experimental data of Ref. [l^ using 
kernel ((67c)) gives — 66 nm, corresponding to 2^ ~ 
11 meV. Analogous analysis was done for the case of 
strongly disordered edge in Sec. IVIIB| and in this model 
one obtains £^ = 120 nm. Indeed, as the disordered edge 
gives more weight to oblique trajectories, as shown in 
Fig. (He), the effective distance from the edge at which 
the kernel decays, is shorter than for the normal inci- 
dence, so a larger value of £^ is required to fit the data. 

The inelastic scattering rate for an electron with the 
energy e due to phonon emission can be written as 
27 = (Ar + A_R-)e/2;^2ji^ where Ar and Xk are dimen- 
sionless electron-phonon coupling constants [Xk is de- 
fined in Eq. (jS.Sp . and Ar is defined analogously, but 
the optical phonons at the T point should be consid- 
ered] . The value of the constant Ar can be reliably taken 
to be about Ar « 0.03. Indeed, a DFT calculation"^? 
gives Ar « 0.028; measurements of the dependence of 
the G-peak frequency ug on the electronic Fermi en- 
ergy ep, duJG/d\eF\ « Ar/27r give Ar « 0.034^1 and 
Ar « 0.027)^ the value of Ar is not renormalized by 
the Coulomb interaction.— The value of Xk has been de- 
bated recently4^i^i^ The measurements of the phonon 
group velocity (see Ref. [5l| for the summary of the ex- 
perimental data) give A^ « 0.04. The ratio between the 
two coupling constants can be also extracted from the 
experimental ratio of the two-phonon peak intensities;^ 
2{XK/Xr)^ w 20,^2 ^hich gives Xk « 0.10. 

Thus, Ar + Xk ~ 0.1 ± 0.03, seems to be a reasonable 
estimate. This estimate gives 27 w 50 meV for electrons 
with the energy uJin/2 — 0.98 eV, which translates into a 
value of £j several times shorter than that following from 
the results of Ref. 



E. On the Tuinstra-Konig relation 

For a sample of a finite size La the total D peak inten- 
sity is proportional to the total length of the edge, i. e., 
the sample perimeter, so that lo oc La- At the same 
time, the intensity of the G peak at 1581 cm^^ is propo- 
tional to the area of the sample, i. e., Iq oc L^. These 
simple facts result in the so-called Tuinstra-Konig rela- 
tion, established in experiments on graphite nanocrystal- 
lites long ago^^*^ and confirmed experimentally many 
times afterwards i^^i^? '^'^ Id/Ig oc 1/La- The propor- 




FIG. 5: Honeycomb lattice with the A and B sublattices and 
the elementarty translation vectors. 

tionality coefficient cannot be determined universally; it 
clearly depends on the character of the edge (as an ex- 
treme case, one can imagine a hexagonal flake with en- 
tirely zigzag edges which do not give any D peak at all, 
except at the junctions between them; then Id is not 
even proportional to La). 

What is the boundary of validity of the Tuinstra-Konig 
relation on the small-size side? At the same time, it 
was noted in Ref. [5^ that since the atomic dispacement 
pattern corresponding to the D peak must involve at least 
one aromatic ring, the size of the ring, a few angstroms, 
represents an absolute lower bound. From the results of 
the present work it follows that the dependence Id o: La 
becomes logarithmically sensitive to the presence of the 
opposite edge, Id oc La ln(wphia/w) if the size is smaller 
the electron inelastic length La < v/{2j), and the whole 
approach becomes invalid when La ~ w/i^ph ~ 4 nm. 
The breakdown of the 1/La dependence has indeed been 
observed for La smaller than a few nanometers;^ 



III. FREE DIRAC ELECTRONS 

In this section we summarize the model for the bulk 
graphene only, which is fully analogous to that of Ref. 4^. 
Properties of the edge are discussed in the next section. 

We measure the single-electron energy e from the Fermi 
level of the undoped (half-filled) graphene. The Fermi 
surface of undoped graphene consists of two points, called 
K and K' . Graphene unit cell contains two atoms, la- 
beled A and B (see Fig. [S]), each of them has one tt- 
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E 


C2 


2C3 


2C6 


0"a,6,c 


/ 

<^a,b,c 




1 


1 


1 


1 


1 


1 




1 


1 


1 


1 


-1 


-1 


B2 


1 


-1 


1 


-1 


1 


-1 


Bi 


1 


-1 


1 


-1 


-1 


1 


El 


2 


-2 


-1 


1 








E2 


2 


2 


-1 


-1 









TABLE I: Irreducible representations of the group C&v and 
their characters. 



irrep 


Ai 


Bi 


A2 


B2 


-El 


E2 


valley-diagonal matrices 


matrix 


1 


A, 




A^Ez 






valley-oflt-diagonal matrices 


matrix 






A^ 


Ay 




-^y'^x ^ -^y^y 



TABLE II: Classification of 4x4 hermitian matrices according 
to irreducible representations (irreps) of the Ce V group. 



results in an additional requirement UtUf = 1. The 
explicit form of Ut depends on the choice of the basis. 

In those rare cases when a specific representation has 
to be chosen, we use that of Ref. |6^ . 

IpAK 

t'^^ , (3.3) 

IpBK' 
_ -IpAK' _ 

where the first subscript labels the sublattice {A, B), and 
the second one labels the valley (K,K'). In this basis 
are the Pauli matrices acting within upper and lower 
2-blocks of the column (the sublattice subspace), while 
are the Pauli matrices acting in the "external" subspace 
of the 2-blocks (the valley subspace). The time reversal 
matrix in this representation is given by Ut = T^yAy. 

The electron Green's function, corresponding to the 
hamiltonian (|3.ip . is given by 



orbital, so there are two electronic states for each point 
of the first Brillouin zone (we disregard the electron spin). 
Thus, there are exactly four electronic states with zero 
energy. An arbitrary linear combination of them is rep- 
resented by a 4-component column vector -0. States with 
low energy are obtained by including a smooth position 
dependence ip{r), r = {x,y). The low-energy hamilto- 
nian has the Dirac form: 



Ho= J (friP^ir) {-ivi:-V)^{r). 



(3.1) 



Here we used the second-quantized notation and intro- 
duced the electronic ?/;-operators ?/;(r), ■0t(r). 

It is convenient to define the 4x4 isospin matrices 
S = (E;r, Ey), not through their explicit form, which de- 
pends on the choice of the basis (specific arrangement 
of the components in the column ip), but through their 
transformation properties. Namely, all 16 generators of 
the SU{4) group, forming the basis in the space of 4 x 4 
hermitian matrices, can be classified according to the ir- 
reducible representations of Cq^, the point group of the 
graphene crystal (Tables III and HI)) . They can be repre- 
sented as products of two mutually commuting algebras 
of Pauli matrices I]a;,Ej,,Ez and A^;, Aj,, A^,^^'^^ which 
fixes their algebraic relations. By definition, are 
the matrices, diagonal in the K, K' subspace, and trans- 
forming according to the Ei representation of Cgy . 

In the following we will take advantage of the sym- 
metry with respect to time reversal. The action of the 
time-reversal operation on the four-component envelope 
function 0(r) is defined as 



7/;(r) ^ Utij*ir), 



(3.2) 



where Ut is a unitary 4x4 matrix, when applied twice, 
the time reversal operation should give an identity, which 



Go(p,e) 



- {vp - i-yeV 



(3.4) 



where p and e are electronic momentum and energy, 
counted from the Dirac point. The inelastic broaden- 
ing 7£ <ti |e| is introduced phenomenologically. In the 
coordinate representation the Green's function is given 

by 



Go(r,e) = 



e -H sgn e — iv 



iI«(C), (3.5) 



where Hq '{(,) is the Hankel function and C = (|e| + 
i^^)r /v. We will mostly need the asymptotic form valid 
at distances r ^ |e|/i', 



Go(r,e) 



sgne + S • r/?' 



1 



8C, 



.(3.6) 



Any wave function ^/'(r) satisfying the Dirac equation, 
(e-t-i7e sgne-f wS- V)?/'(r) — 0, in some region O of space, 
satisfies also the Huygens-Fresnel principle. Namely, the 
value of iP{y) at an arbitrary point r G O can be written 
as an integral over the boundary 90, 

V'(r) = iv j n-[Go{Y~ r^, e) S ^j{r,)] dr,. (3.7) 

do 



Here is the distance along the boundary, n is the inner 
normal to the boundary. This relation follows from the 
Gauss theorem and the fact that {p^ + V'^)HQ^\pr) = 
4i(5(r) for any p. 
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IV. MODELS FOR ELECTRONS NEAR THE 
EDGE 

A. Translationally invariant edge 

The main assumption of this subsection is that the 
component of the electronic momentum p along the edge 
is conserved upon reflection, so that a plane wave is re- 
flected as a plane wave. The most studied ideal zigzag 
and armchair edges fall into this category. Here we do 
not restrict ourselves just to zigzag or armchair edges, 
requiring only that the spatial period de of the edge is 
smaller than half the electron wavelength, de < TrAg. For 
de ^ the reflection of electrons from the edge can 
be described by an effective boundary condition for the 
electronic wave functioni^i^l 

The edge is assumed to be a straight line determined 
by its normal unit vector n, so that graphene occupies 
the half-plane n • r > 0. The microscopic Schrodinger 
equation determines the effective boundary condition on 
the wave function ijjir), which for smooth functions (on 
the scale de) can be simply written as Bijj\^^^^ = 0, 
where is a 4 x 4 hermitian matrix. The rank of B 
is equal to 2 since the linear space of incident states 
at fixed energy is two-dimensional due to the valley de- 
generacy. Thus, it has two zero eigenvalues, while the 
other two can be set to 1 without the loss of generality 
(only the zero subspace of B matters). Thus, one can 
impose the condition B^ = B. Equivalently, one can 
write B = (1 — M)/2, where M has the same eigen- 
vectors as B, but its eigenvalues are equal to ±1, hence 
= 1. To ensure current conservation, the condition 
BV'lodgo ~ ^ niust automatically yield ijj^n ■ S)?/; = 0; 
this means that Af(n • E) + (n • S)Af = 0. Finally, 
the time reversal symmetry requires that the conditions 
Bip — and BUt^p* — must be equivalent, which yields 
M* = u} AlUt- To summarize all the above arguments, 
the general energy-independent boundary condition has 
the form 

(1 - Am^^^^ = 0, (4.1) 

where the 4x4 hermitian matrix M satisfies the following 
conditions, which result to be the same as obtained in 
Ref. M- 

= 1, M(n • S) -t- (n • S)A/ = 0, Af = UtM*uj . 

(4.2) 

Matrices satisfying these constraints can be parametrized 
by an angle x ^^'^ ^ three-dimensional unit vector 
{nix, my, m^)-^ 

M = (S^ cosx + [n X S]^ sinx) A//a, (4.3a) 
A^A = ^ m^A„ m, e R, ^ mj ^ ];4.3b) 

i—x,y,z i—x,y,z 

Without loss of generality we can assume cos x > (the 
negative sign can always be incorporated in M\). 



Explicit expressions for the matrix M, corresponding 
to the edges shown in Fig. [TJ can be obtained in the tight- 
binding model on the terminated honeycomb lattice^'^ [it 
is convenient to use representation (13. 3p ]. For the zigzag 
edge with n = —ey [Fig. [TJa)] the boundary condition is 
i'BK = ^BK' = 0, which gives M — SzA^. This agrees 
with the prediction of Ref. [l^ that upon reflection from 
a zigzag edge the electron cannot scatter between the 
valleys, which corresponds to a valley-diagonal matrix 
M. For the armchair edge with n = e^, [Fig. [T](a) rotated 
by 27r/3 counterclockwise] we have tpAK + ipAK' = '4'bk + 
i^BK' — 0, and M — ^T,yAy. It can be shown that in 
the nearest-neigbor tight-binding model on a terminated 
honeycomb lattice only zigzag and armchair boundary 
conditions can be obtained, the latter occurring only if 
the edge direction is armchair, while to obtain the full 
form of Eq. (|4.3a|) , one has to include an on-site potential 
in the vicinity of the edge.^^ 

It is known for quite some time that a perfect zigzag 
edge supports states, confined to the edgei2ii22ii^i2ii2^ Let 
us see what class of boundary conditions is compatible 
with existence of such states. The wave function of an 
edge state must have the form 

^(r) = Voe"''^""")+'''"["'""l% ^>o, (4.4) 

where the vector V'o is such that the solution satisfies both 
the Dirac equation in the bulk with some energy e, as well 
as the boundary condition at the edge. It is convenient 
to make a unitary substitution ipo = e^^'''^/^e~^^''^"/^^o, 
where ifn is the polar angle of the direction n. Then the 
two conditions have the following form: 

(ikS^ -I- p|| cos x'^y + P\\ sin x^z) ^ ^ V'o(4.5a) 

S.A^A^o = ^0- (4.5b) 

The boundary condition is satisfied by two linearly in- 
dependent vectors V± which can be chosen to satisfy 
^zV'i = iV'ib (to find them it is sufficient to diagonalize 
the matrix AfA). Each of them satisfies the first condi- 
tion if and only if e = ±i;p|| sinx, k = =Fp|| cosx- The 
requirement k > leaves only one of them, ■(/'_ sgnpn 7 and 
the energy of the edge state is e = | sinx- Thus, it 

seems that almost any edge can support a bound state, 
the exception being the case cosx = (armchair- like 
edge), which thus seems to be a special rather than a 
general case. 

Now we turn to the scattering (reflecting) states, which 
are the ones responsible for the edge-assisted Raman 
scattering. Even though the component py of the electron 
momentum p, parallel to the edge, is conserved, reflec- 
tion can change the valley structure of the electron wave 
function. The general form of such solution is 

l/,(r) = ^pe'P^(n-r)+»P||[nxr], _^ 

+ S'A[n X S];,i/;pe-^P^("-'")+*fi|["^'"l% (4.6) 

where p± — {n ■ p) < 0, = [n x p]^, and Vp is an 
eigenvector of (p • T^)ipp — ±\p\'iljp. The first term rep- 
resents the wave incident on the edge, the second one is 
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the reflected wave. The matrix [n x SJ^ simply aligns 
the isospin of the reflected particle with the new direc- 
tion of momentum. The unitary matrix Sa represents a 
rotation in the valley subspace. It should be found from 
the boundary condition (j4.ip (this is conveniently done 
in the basis of the eigenvectors of A/a), which gives 



Ao + Ma(cosx + C* sinx) 
Ao + Ma (cos X + C sin x) ' 



( = ±^{-[n X p], + i{n ■ p)} , 



(4.7a) 
Id = 1. (4.7b) 



For sin X = (zigzag edge) we have Sa — (■ For cos x = 
(armchair edge) we have S'a — Ma, independent of the 
direction of p. The reflected part of Eq. (|4.6p can be 
identically rewritten using the Huygens-Fresnel principle, 
Eq. dsn, as 

V'(r) = Vpe'P"" - J dreGo(r-re,e)«S,5AVpe'P'"% 

edge 

(4.8) 

so that —v'SzSa can be viewed as the T-matix of the 
edge. 

When S'a does not depend on the direction of p (arm- 
chair edge), it is easy to write down the exact explicit 
expression for the single-particle Green's function: 

G(r,r';e) = Go(r-r',e) 

+ Go(r - r' + 2n(n • r'), e) [n x SJ^Sa. 

(4.9) 

The second term represents nothing but the contribution 
of a fictitious image source of particles, appropriately ro- 
tated, and placed at the point r' — 2n(n • r') obtained 
from r' by the reflection with respect to the edge. In 
the quasiclassical approximation (analogous to geometric 
optics), Eq. (|4.9p is also valid for a general edge at large 
distances r,r' 3> v/\e\, provided that position-dependent 
S'a is taken, determined by 



C = - sgn e 



z(n-(r + rO) + [nx (r - rQ]^ 
v/(n.(r + r'))2 + [nx (r-r')]r 



(4.10) 



Again, using the Huygens-Fresnel principle, we can 
rewrite Eq. (14. 9|) identically as 



G(r,r';e) = Go(r-r',e) 

J dYe Go(r - re, e) wS^Sa Go(re - r', e) 



edge 



B. Atomically rough edge 



(4.11) 



As discussed in Sec. Ill A[ when the edge is rough on the 
atomic length scale, electron scattering is random both 
in all directions and between the two valleys. This case 



will be of main interest for us, as it (i) represents the 
opposite limiting case to that of an ordered edge, and 
(ii) can be described by a simple model proposed below. 
The main assumption is that each point of an atomi- 
cally rough edge acts as an independent point scatterer, 
independent from other segments. Associating thus a T- 
matrix to each point of the edge, we write the scattered 
wave function in the form 

^(r) = Vpe'P"- 

+ j dreGo(r-re,e)r(s,ep;re)Vpe'P(;4,12) 

edge 

'-|r-re|' ^--|p|' 

where energy argument of the Green's function e = 
±up for electrons and holes, respectively. The (one- 
dimensional) integration is performed along the edge, 
which is assumed to be a straight line determined, as 
in the previous subsection, by the condition (n • rg) = 0, 
where the unit vector n is the normal to the edge. The 
unit vectors Gp and s indicate the incident and scattering 
directions. 

The T-matrix must satisfy (i) the particle conservation 
condition (unitarity) , and (ii) the time reversal symmetry 
(reciprocity). 



r(s, ep; re) = Ut T^{~ep, -s; re) Ul 



(4.13) 



Here stands for the 4x4 matrix transpose. Unitar- 
ity and reciprocity are discussed in Appendix \^ in the 
context of a general scattering theory, similar to that for 
light scattering on a rough surface?2i We propose the fol- 
lowing form of the T-matrix: 



r(ps,p;re 



p(s,ep)t;E3SA(re). 



(4.14a) 



The angular factor p(s,ep) ensures the particle conser- 
vation (unitarity of the edge scattering, see Appendix \X\ 
for details). 



P(s,ep) 



-2(n • ep)(n • s) 



1 - (n • ep)(n • s) - [n x ep]z[n x s]^ ' 

(4.14b) 

If we introduce the angles of incidence and scattering by 
(n • Bp) = - cos (pi, (n • s) = cos(ps, [n x ep]^ = sirnpi, 
[n X s]^ — simps, the specular direction corresponding 
to ips = ifii, then p — 2 cos cos (^s/[l + cos{(pi + ips)]- 
Note that since the structure of the wave functions in 
the S-subspace is fixed by the direction of momentum, 
one may suggest slightly different forms of Eq. (j4.14bp 
which would be equivalent. For example, the T-matrix 
obtained by from Eq. (j4.14bp by replacement of by 
i[n X S]z and by changing the sign of the third term 
in the denominator in Eq. (I4.14bp would have the same 
matrix elements between Dirac eigenstates of the same 
energy. 
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For a coordinate-independent the re integration 
eliminates all directions s different from the specular 
one. Since for the latter p = 1, Eq. (|4.14ap reduces to 
Eq. (|4.8p . In this case can be identified with the scat- 
tering matrix. For the short-range disorder such identi- 
fication is not possible, since scattering process is nec- 
essarily non-local on the length scale of at least ^ 1/p. 
The connection between the T-matrix in Eq. (I4.12p and 
the scattering matrix is more complicated, and is discuss 
in detail in Appendix |^ Here we only mention that it 
would make absolutely no sense to require the unitar- 
ity of S\{re) at each given point. Instead, we use the 
following form: 



S'A(re) = n7(re) ^ mj(re)Ai 



m,2(re) = 1. 



E 

i—x.y,z i—x,y,z 

(4.14c) 

Here in (re) is a complex gaussian random variable whose 
real and imaginary parts are distributed independently 
and identically (so its phase is uniformly distributed be- 
tween and 27r). The numbers ma:,my,mz, which can 
be viewed as components of a unit three-dimensional vec- 
tor, must be real to ensure the time-reversal symmetry. 
One may assume them to be constant or taking just few 
definite values, which would correspond the edge to be 
composed of segments with a <C de ^ 1/p of definite 
types (e. g., zigzag or armchair). For de a, when the 
scattering between the valleys is completely random, the 
vector (mx,my,mz) can be taken uniformly distributed 
over the unit sphere. We assume the matrices S'A(re) to 
be uncorrelated at different points on the edge, distant 
by more than de, by writing 



wire) w*(r'e) = -g- S{re - r'J. 



(4.14d) 



Here the overline denotes the ensemble averaging. How- 
ever, we assume that this product is self-averaging upon 
spatial integration, i. e., that Eq. (j4.14dp holds even in 
the absence of the ensemble averaging when integrated 
over a sufficiently long segment of the edge (namely, 
longer than de). The prefactor 7ri>/|e| at the 5-function 
ensures the unitarity of scattering, see Appendix [XI 

Eq. (|4.12p for the wave function yields an analogous 
expression for the Green's function, valid sufficiently far 
from the edge, (n • r) ^ '^/kli (n ' r') ^ v/\e\: 

G(r,r';e) -Go(r-r',e) + 

J dve Go(r - re, e) T(s, s'; re) Go(re - M) 



edge 

r — rp 



_ <-e I 

^ ~ I I' ® ^ Iv' 



r re 



(4.16) 



V. PHONONS AND RAMAN SCATTERING 

We restrict our attention to scalar phonons with wave 
vectors close to K and K' points - those responsible for 



? 



r 



r 



A, B, 

FIG. 6: Plionon modes responsible for the Raman D peak. 



the Raman D peak. The two real linear combinations of 
the modes at K and K' points transform according to Ai 
and Bi representations of Cqv and are shown in Fig. [B] 
We take the magnitude of the carbon atom displacement 
as the normal coordinate for each mode, denoted by Ua 
and Ub, respectively. Upon quantization of the phonon 
field, the displacement operators Ua,Ub and the lattice 
hamiltonian 7?ph are expressed in terms of the phonon 

a, 6, as 



creation and annihilation operators b„ 



Hph — LxLy 



d\ foq^e'i-- 



(2^)2 ^2NMijph 
rf2q ^ 
(2^ "P'^ 



(5.1a) 



1 (^^q^^q/^ + \ 



15. lb) 



The crystal is assumed to have the area L^Ly, and to 
contain N carbon atoms of mass M . The area per carbon 
atom is L^Ly/N = y/27a'^/4:. 

The phonon frequency ojph ~ 1350 cm~^, standing in 
Eq. (|5.1bp . is assumed to be independent of the phonon 
momentum. To check the validity of this assumption 
one should compare the corresponding energy scale (the 
spread of the phonon momenta Ag multiplied by the 
phonon group velocity Wph) with the electronic energy 
uncertainty. The latter is given by Wph itself. Recall- 
ing that phonon emission corresponds to the backscat- 
tering of the electron (hole), Aq is given by the uncer- 
tainty of the electronic momentum, Aq ujph/v. Since 
Vph/v « 7-10^'^ ^ 1, the phonon dispersion can be safely 
neglected. 

If we neglect the phonon dispersion, the normal modes 
and the phonon hamiltonian can be rewritten in the co- 
ordinate representation by introducing the creation and 
annihilation operators for a phonon in a given point of 
the sample: 



LxLy 



/ 2NMujpi, 



4^(r) + $t (r) 
<&M(r)<&J.(r) + 



N 



2LxLy 



(5.2a) 

(5.2b) 
- ;5.2c) 
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Then it is convenient to define the phonon Green's func- 
tion as the time-ordered average of the ^-operators, 



4+)(r,c.) = I (T<|.^(r,t)<|t(O,0))e'"*d2rdt 
Sir) 



u — Wph + io 



(5.3) 



By symmetry, in the electron-phonon interaction hamil- 
tonian the normal mode displacements couple to the 
corresponding valley-off-diagonal 4x4 matrices from Ta- 
ble HIl^^ 



K 



d'ri^^r) K(r)A,S, + Ufc(r)Aj,E,] i^(r). 

(5.4) 

Here Fk is the coupling constant having the dimension- 
ality of a force. It is more convenient to use the dimen- 
sionless coupling constant 



A 



K 



P2 



(5.5) 



The value of Xk was discussed in Sec. IIIDI 

The hamiltonian describing interaction of electrons 
with light is obtained from the Dirac hamiltonian (|3.1[) 
by replacement V — > V — z(e/c)A, where the vector po- 
tential A is expressed in terms of creation and annihi- 
lation operators 0'Qg,o,Q,e of three-dimensional photons 
in the quantization volume V — L^LyLz, labeled by the 
three-dimensional wave vector Q and two transverse po- 
larizations ^ = 1,2 with unit vectors eq^^: 



A(r) = E 



27rc 
VQ 



iQr 



h.C. 



(5.6) 



The derivation of the formal expression for the Raman 
scattering probability is fully analogous to that given in 
Ref. m. The only difference is that the calculation is done 
in the coordinate representation. As a result, we obtain 
the following expression for the probability for an inci- 
dent photon with wave vector Qi„ and polarization ei„ 
to be scattered with emission of a single phonon within 
an elementary area d^R around a given point R: 



dip 



CLxLy 



27r(5(c|Q, 



^out) 

(5.7a) 



2 ^l-OinUloiit 



2^ 



X e 



'TrjP^-l-P^}, (5.7b) 
G( ^){Gin ■ S) G(ri„, R; e — uiin) 

X A^^yT,^ G{R, rouu e - ^out){e*o^t ' S),(5.7c) 
G( e) (eL,.S)G( 

I out I R; e + ^out) 



X K^^yY.^ G(R, Yin]t + UJin){ei, 



(5.7d) 



Here G(r, r';e) is the electronic Green's function corre- 
sponding to the full single-particle part of the hamilto- 
nian (i. e., including not only the Dirac term, but the 
edge as well) . It can be represented in terms of the exact 
single-electron eigenfunctions t/'s (r) and energies as a 
sum over the eigenstates s: 



G(r,r';e)=E-^ 



V'.(r)^Kr) 



i^e sgn e 



(5.8) 



Using this representation, integrating over the energy and 
the coordinates one obtains Eq. (|2.ip . The summation 
in Eq. I|5.7ap is performed over the wave vectors Qout 
and the polarizations Gout of the scatered photon. When 
integrated over the area of the crystal, Eq. (|5.7ap gives 
the absolute dimensionless probability of the one-phonon 
Raman scattering for a single incident photon. 

The matrix element M^^ can always be represented in 
the form 



>^M=>'u(eLt)*+M;i(eLi)* 



(5.9) 



If one collects all the light scattered in the full solid angle 
47r, without analyzing the polarization, the integration 
over the angles of Qout is straightforward. It gives 



dl 



D 



d2R 



E i 

^ fj.=a,b 



\Ml\') 



(5.10a) 



The dependence on the polarization of the scattered light 
is obtained most easily when the light is collected in a 
small solid angle Oout 47r around the normal (the case 
of an arbitrary solid angle was considered in Ref. H^)- If 
the analyzer is oriented at an angle ipout to the x axis, 
the polarization-dependent intensity is given by 



dip 
d^R 



Oout 



C LxLy _ 



LMJj cos if out + Ml simpo 



a.b 



(5.10b) 

Eq. (|5.6p corresponds to the free-space quantization of 
the electromagnetic field whose normal modes are plane 
waves. In the case of a spatially resolved experiment as 
in Ref. Tg", Eqs. (|5.7ap . (j5.7bp in order to account for the 
spatial profile of the electric field, induced by the focusing 
lens. Namely, the electric field, corresponding to a single 
photon with a wave vector Qm, incident from vacuum, 
should be replaced by the field 5m(i"m) of the focused 
laser beam. 



/ 27ra;i„ 

LxLyLz 



(5.11) 



As long as the distance between the lens and the sample 
is much larger than the light wavelegth, the summation 
over the continuum of the final states of the scattered 
photon can still be performed using the vacuum mode 
structure. Finally, dividing the resulting probability by 
the photon attempt period L^/c, we obtain the number 
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of photons emitted per unit time, dlu /dt (which is more 
appropriate when the incident hght is characterized by 
its electric field strength). As a result, Eqs. (|5.7ap . (|5.7bp 
are modified as follows: 



dip 
d^Rdt 



1 



27rwf 



E 



Gout 



out 

(27r)3 



27r(5(c|Q„ 



X \2M 

fi—a,b 



2 

Ml ' 



■ l-Oout) 

(5.12a) 



de 



■d'^r d'^r 



2 ^'-^in^out J 27r 

X £„(r„)e-*'2°-'^°-Tr{l?^ +P,,](5.12b) 



VI. RAMAN SCATTERING ON A 
TRANSLATIONALLY INVARIANT EDGE 

A. General considerations 



For simplicity we consider an armchair edge charac- 
terized by n — e^, M — —Y^yKy, as discussed in the 
previous section. Due to the translational invariance in 
the y direction, it is sufficient to calculate the probabil- 
ity of phonon emission at the point R = Xe^.. As we 
will see below, the main contribution to the signal comes 
from X 3> Ae = 2v/ujin. In this regime the motion of the 
photoexcited electron-hole pair can be described quasi- 
classically, and the asymptotic large-distance expansion 
for the Green's functions, Eq. p.6p . can be used. Namely, 
the electron and the hole can be viewed as wave packets 
of the size ~ ^/X^X, propagating across the crystal along 
classical trajectories. Initially, they are created in the 
same region of space of the size ~ \/XeX with opposite 
momenta and opposite velocities. As they undergo scat- 
tering processes (emission of a phonon or refiection from 
the edge), they change the directions of their momenta. 
In order to recombine radiatively and contribute to Ra- 
man signal, the electron and the hole should meet again 
within a spatial region of the size ~ ^/X^X. Clearly, these 
conditions can be fulfilled if all momenta are almost per- 
pendicular to the edge. Small deviations by an angle 
~ \/X^JX are allowed by quantum uncertainty. These 
considerations are illustrated by Fig. [3l 

Emission of one of the phonons shown in Fig. [S] cor- 
responds to intervalley scattering of the photoexcited 
electron or the hole, as represented formally by one of 
the valley-off-diagonal matrices A^; or Aj^ in the ma- 
trix element, see Eqs. (j5.7bp - (j5.7dp . For the process to 
be allowed, another act of intervalley scattering is re- 
quired, so one of the three electronic Green's functions 
in Eqs. (|5.7cp . (j5.7dp must contain another A^; or Aj, ma- 
trix from the decomposition (|4.9p . Otherwise, the trace 
in Eq. (j5.7bp vanishes. Thus, for an armchair edge with 
M = —KyTiy only the Bi phonon can be emitted, so 
M.a = 0, Mb — M.. Trajectories, corresponding to 
decomposition of each of the three Green's functions in 



Eq. (|5.7cp are shown in Fig. [3] (a), (b), (c), respectively. 
According to the quasiclassical picture, the electron and 
the hole have to travel the same distance between cre- 
ation and annihilation, as their velocities are equal. Then 
the process in Fig. [3] (a) has more phase space satisfying 
this restriction, and this process gives the main contri- 
bution to the Raman matrix element. This will be also 
shown explicitly below. 

The general polarization structure of the matrix ele- 
ment, compatible with the reflection symmetry y — > —y 
possessed by the armchair edge, is 



M 



M\\ Sirnpin SinLpout + Ma_ COSLpm cosyjoUS-l) 



(we introduced (pmjfout, the polar angles of the polar- 
ization vectors). Since the interband transition dipole 
moment is perpendicular to the electronic momentum, 
and for a regular edge the momentum is (almost) per- 
pendincular to the edge, A4± is entirely due to quantum 
diffraction. Thus, A4± is smaller than A^|| by the pa- 
rameter Xc/X <^ 1 (it is this parameter that governs the 
quantum diffraction, as discussed above). Nevertheless, 
in a polarization-resolved experiment the two intensities 
can be measured independently, so below we calculate 
both A^ii and Mj_, each to the leading order in X^/X. 



B. Spatially integrated intensity and polarization 
dependence 



In this subsection the electric field of the excitation 
wave is assumed to be spatially homogeneous. As dis- 
placements along the edge are expected to be parametri- 
cally smaller than those away from the edge, \y\ ^ a/JAT, 
we use the paraxial approximation, |r| « |a;| -I- y'^ /{2\x\). 
The Green's function can be approximated as 



Go(r,e) 



^{i\e\--t)\x\/v+i\e\y^/{2v\x\) 



27rt;'^|a::| 

sgn e -I- Sj, sgn x 



2\x\ 



(S- 


iv ^ 


1 




\ex\j 





(6.2) 



where the coefficient at each matrix is taken in the lead- 
ing order. Taking the first term in the square brackets 
in Eq. (|6.2p and evaluating the matrix traces, we obtain 
the following expression for A^n, corresponding to the 
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process in Fig. [3] (a): 



Mn 



oo X 



j ^ j dx^ndXout j dyindyoute''^ 2iX/v 

(6.3a) 



-oo 



[ujin - t){LOout - e)e| 



V 2v{xin + XoMt) 



V 2v{X~Xin)' 



■ .3b) 



For the calculation oi A4± we need the rest of the terms 
in the square brackets in Eq. (|6.2p . As a result, TWx is 
given by an analogous integral, but with an extra factor 
in the integrand. 



yout Uin 
•^in ~t~ Xout 

iv 



X 



IJout 



Vout 



^out 



4e U 



•^out ) 



The details of integration are given in Appendix|B] First, 
we integrate over y^ and yout- The subsequent integra- 
tion over e fixes Xm/v + Xout/v ~ (X — Xin)lv + {X — 
Xout)/v [the difference is allowed to be ^ \/X/ (uwi„) <C 
X], which has the meaning of the time spent by the 
electron and the hole in traveling from the creation 
point Xin to the annihilation point Xout- At the same 
time, x-integration fixes e « uJmj^ with the precision 
~ \J uJinV / X <C oJin- Performing all the integrations, we 
obtain the matrix element, 



M 



ttXk 



[u;phX/{2v)] 



X e 
iM 



iuji„X/v Wph/ (2w) 

i{uJi„+uJa^t)X/{2v)-2'yX/v 



(6.4a) 
(6.4b) 



According to Eq. (|5.10ap . the integrated intensity into 
the full solid angle 47r, summed over two polarizations of 



the emitted photon, is given by 



d2R 



8Xk (eVc)2 col sin2[c.phX/(2«)] g-^-r^/" 



3 LxLy 
X [sin^ ip^n - 



c2 Kh/(2z;)]2 



UinX/v 

(6.5a) 



8A 



K 



ID 



sin^ ipi. 



In 



(47)' 



COS^ ifi, 



(47)2 



-2-ln- 



t^ph 



(6.5b) 



The second term in the square brackets is written with 
logarithmic precision, since the short-distance cutoff ^ 
v/ujin is known only up to a factor of the order of 1. If 
we use Eq. (I5.10bp for the intensity in a solid angle Oout 
in the presence of an analyzer, we obtain 



4A 



K 



Oout 

An 



^ph 



sin^ ifiin s\T? ifout In 



^ph 



(47)2 



-I- COS^ (/Jj„ COS^ ipout 



(47)2 



'^ph 



.(6.5c) 



Let us estimate the contribution to the matrix ele- 
ment A4' corresponding to Fig. [3](b), i. e., when decom- 
position (|4.9p is applied to G(R, ro„t; e — t^out)- Eq. (|6.1| 
remains valid, as it is based on symmetry only. The ex- 
pression for A^Ji looks analogous to Eqs. ()6.3ap . (|6.3bp 



M[ = — 



iX 



K 



de 




dXn 



\{u!in - e){u!out - e)e| 



dy^ndyout e'*'"^^^/- x 
(6.6a) 



iX-X^n){X + 

^out ){Xrn - Xout) 
\e\{Xin - Xout) , \e\iyin - yout)'^ , 



V 

\0Jout ■ 



2v{xi 
e\{X + Xout) 



Xout ) 
\l^out 



^IvLt 



\uJin ~ e\iX - Xtn) , I 



2v{X + Xout) 
,2 



2v{X - Xin) ' 



(6.6b) 



However, here integration over e fixes Xm —Xout ~ Xout + 
X + (X — Xin), which is compatible with the limits of the 
spatial integration only when Xout ^ v/ujin, X ~ Xin ^ 
v/bJim &s shown in Fig [3] (b). This restriction results in 
suppression M',JM\\ ^ uu^^h/uin 1. 
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If Xout, \X — Xin\ ~ v/ujin, the asymptotic form, 
Eq. (|6.2p . cannot be used for G'(rj„, R; aj,„ — e) [thus, 
Eqs. (|6.6ap . (|6.6bp represent only an estimate of A^J| 
by the order of magnitude], but it can be used for 
G(R, Touui^out - e), and G{rout,rin; e). This fact results 
in an additional smallness for the matrix element A4'^: 
assuming the typical value X ~ u/wph, we can write 
lA^'^l - |M|||(t^phM„) - |A^i|(L^phM„). Thus, M'^ 

produces only a small correction to the cos^ (pin term in 
Eqs. (I6.5a|) . (I6.5b|) . and to the cos^ (fin cos^ ^pout term in 
Eq. ((63cl) . 

Finally, the intensity in Eq. (|6.5cp has an interfer- 
ence contribution cx Re{A4J^A4'^), which produces a term 
oc sin ipin cos (fin sin ipout cos ifiout ■ Note that the interfer- 
ence term Re(A^jjA^j_) is absent because of the factor i 
in Eq. (|6.4b[) . We have not been able to calculate A4'^ 
explicitly or to establish a phase relationsip between 
and Ai'j^ in any other way. However, it is hard to imag- 
ine that the inteference of two strongly oscillating am- 
plitudes, corresponding to two different processes, would 
survive the integration over X. 



C. Excitation position dependence 

This subsection aims at describing a spatially resolved 
experiment like that of Ref. [l^ and clarifying the role of 
different length scales in the dependence of the Raman 
intensity on the position of the excitation spot. Con- 
sequently, we use Eqs. (I5.12ap . (j5.12bp . and repeat the 
calculation of Sec. IVI Bl with an arbitrary dependence of 
fi„(r), smooth on the length scale v/ujin (we assume de- 
tection without analyzer and sum over the polarizations 
of the scattered photon). The result is 



-rr [ — ] - sin ifiy 

at 3tt \ c J c 



d2/X(— ,46.)7fi) 

V^ph 27/ 



I{eph,£y)^ J £in{x,y) £*n{x' ,y) IC{x,x') dx d3(^.7h) 


K:ix,x') = -g-^^^-^'^/^p" Ei(-2max{j;,j;'}/^^)(6.7c) 
where the exponential integral Ei(z) is defined as 



00 

Ei(-z) = J 



t 



(6.7d) 



excitation profile w{x) is assumed to be much longer than 
^ph = v/ujpii and the electron inelastic scattering length 
£~i — v/{2j): £ph,i-y ^ L. In all above expressions of 
this section no assumption was made about the relative 
magnitude of iph and £-^. However, it is reasonable to 
assume iph ^ if] also, the final expressions become more 
compact in this limit. 

In the zeroth approximation in l/L one can assume 
w{x) = const in the effective region of integration in 
Eq. (j6.7bp . i. e. replace the kernel by a (5-function. This 
gives the result of Sec. IVIBl 



icco^iM-^o) + o{e/L), 



Iq — IC{x,x') dx dx' = ipi^ln 



+ 4^2 



(6.8a) 
^ (6.8b) 



The length scale Zq, appearing here, may be viewed as the 
effective range of integration in Eq. (|6.7bp . which deter- 
mines the magnitude of the signal. As we see, this length 
is mainly determined by i'ph and is only logarithmically 
sensitive to the electron inelastic scattering. 

What is detected experimentally in Ref. [la is the differ- 
ence in the profiles of u)(— xq) and X^g , appearing because 
of the non-locality of the kernel. This difference appears 
in the second order of the expansion of Eq. (feJbl) in the 
spatial derivatives of £in{x) (i. e, in the order l/L^): 

I.o = - xo) + ^ w"{h - xo) + 0{e/L^). (6.9) 

Here the length [the "center of mass" of the kernel in 
Eq. (|6.7cp ] is given by 



Re / xlC(x,x') — ^2 — 



(6.10) 



It describes the overall shift of the profile, which may be 
difficult to detect experimentally, unless the precise loca- 
tion of the edge is known. The length /2, appearing in 
Eq. (|6.9p , determines the broadening of the signal profile 
Xj;„ with respect to the excitation profile w(— xq), pro- 
portional to w" (x) (the second derivative). In the limit 
i-y ^ £ph it is given by (see Appendix [Cl for the full ex- 
pression and other details): 



Let us assume the excitation intensity to have the form 
|^^m(a;)|2 = w{x — xq), where w{x) is some smooth bell- 
shaped function centered at a; = 0, and xq is the position 
of the focus, which serves as the experimental control pa- 
rameter. In the following we also assume that the phase 
of £in{x) does not depend on x, then £in{x) can be taken 
real without loss of generality. The integral in Eq. (|6.7bp 
is determined by three length scales. The width L of the 



;2^^2 21n(^^/^ph) 



161n'(^^/^ph) 



0(^^j_ 



(6.11) 



Note that this length is indeed determined by the elec- 
tronic inelastic length (up to logarithmic corrections), in 
qualitative agreement with the assumption of Ref. jig . 

Instead of Eq. (|6.7cp . Cangado, Beams and Novotny^^ 
fitted the experimental profile T^o using the following ex- 
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pression: 



(6.12) 



where the excitation profile was independently deter- 
mined to be eaussian: 

\£,n{x)\'^ oc 6-=^'/-^'. The effec- 
tive position of the edge Xe, the width as well as the 
overall proportionality coefficient were used as fitting pa- 
rameters, and the value = 20 nm was obtained. Ex- 
panding Eq. (|6.12|) to the order 1/L^ and comparing it 
to Eq. (IC5ap . we obtain 



Ae- 



0{L-'). 
(6.13) 



This equation is satisfied for all Xq provided that Xe 
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Thus, in spite of 



XD = h, A ^ I + Ij/L"^, 

the fact that in Ref. [l^ a wrong kernel was used, we 
still can take the experimentally found I2 and use it with 
the correct kernel. Namely, the experimentally measured 
xd — 20 nm yields I2 — 14 nm. Using Eq. (|6.1ip and 
taking = 4 nm, we obtain £j = 66 nm, which gives 
27 = 11 meV. As discussed in Sec. HID I and in Ref. iH, 
this value is significantly smaller than an estimate ob- 
tained using other sources of information. 



VII. RAMAN SCATTERING ON AN 
ATOMICALLY ROUGH EDGE 



In this section we calculate the Raman intensity for an 
edge rough on atomic scale, and described by the model 
of Sec. IIVBI The general arguments of Sec lVIAl mostly 
remain valid, except for that of the symmetry y —* —y, 
not possessed by any given realization of the disorder. 
This symmetry is restored upon averaging of the inten- 
sity over the realizations of disorder, but the matrix ele- 
ment Mfj^ must be taken in the general form. 




FIG. 7: (Color on-line.) Electron trajectories corresponding 
to Raman scattering on a disordered edge. Notations are the 
same as in Fig. [3] 



ing expression for the Raman matrix element: 




X / dy,e-2(^/'^)V^^Tr{A^5A(ye)} X 



(X2 + 2;2)3/2 

d'£,ind^out / dijindr/out X 




\ujin - e\\uJout - ellep 



'S (7.1a) 



2 
Va 



C.-j,(7.1b) 
(7.1c) 



where a is either "in" or ''out" . Analogously to the case 
of the regular edge, we first integrate over rjin,r]out, and 
subsequent integration over e and ^in,out fixes ^in+^out ~ 
^A2 + yl^ € « uJin/2. As a result, we obtain 



Spatially integrated intensity and polarization ^g2 j 

dependence Mf, = "T V ^ / ^^^e Tr {A^5A(ye)} X 



Since the edge can scatter an electron by an arbitrary ^ ^(iuin/2+iu}o^t/2-2-^)^Jx^+yl/v ^ 
angle, it is convenient to use the rotated coordinates YiY^cy s^n \(Ycy _i_ 

(4,77), as shown m lig.LlI lakmg the first Green s func- 
tion in Eqs. ([g?7c)l . (|5.7dp as given by Eq. (jiTTC)) . and 



(A2 + y2)2 



taking the free Green's functions in the paraxial approx- sin[a;ph\/ A2 -|- y'^/(2v)] 

imation (i. e., assuming jr/j <C £,), we arrive at the follow- ujp\^/{2v) ' \ ■ ) 
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To calculate the intensity, we use Eqs. (I4.14cp . (j4.14dp 
to average the square of the matrix element, and sum 
over the two phonon modes. Angular integration and 
summation over the two detector polarizations according 
to Eq. (|5.10a|) gives 



dip 



2 ^ 



sin^[c^phVX2 + j/g/(2i;)] 
X^{X^ sin^ fin + yl cos^ tpm) 



36 



(47)^ 



(47)2 



X (Ssin^ ipin + cos^ Lpin) 



(7.3a) 



(7.3b) 



Eq. (|5.10b[) for the intensity emitted into a solid angle 
Oout iu the presence of an analyzer gives 



In = 



24 47r 



c2 wL^ ujI^ (47)2 
1 



(47)2 



sin + sin ipout + ^ cos(2(^i„ - 2(po„4) 

(7.3c) 

The trigonometric expression in the square brackets 
can be identically rewritten as sin2 Lpi^ + sin(2(y5o„t — 
LPin) sin(^i„ + l/2, so its absolute minimum is 1/4, reached 
at i^out = -<^m = ±7r/6. 



B. Excitation position dependence 



Here we follow the same logic as in Sec. IVI C[ but in- 
stead of Eqs. (|6.7ap - (|6.7cp we have 

dip _ / 
dt dir \ c I c 



X 



X^{X sin (p^n + ye cos iy9i„)2 





c I £.inX Cm2/e 
X tin I ,Y + ye 



where (X, F) is the point where the phonon is emitted. 
As in Sec. IVI CI we expand in the spatial derivatives of 
£in{r), and obtain Eq. (|6.9p with I2 given by (actually, 
the result depends slightly on the polarization; here we 
choose unpolarized detection and excitation polarization 
along the edge, (fin — ■7t/2): 



51n(£-^/£ph)-3-2^V(457r) 
' 481n2(£^/^ph) 



0(^2 



(7.5) 



This expression reproduces the experimental value I2 
14 nm if = 120 nm (again, £p\i = 4 nm is taken). 



VIII. RAMAN SCATTERING ON A 
FRAGMENTED EDGE 



In this section we consider the Raman scattering on 
an edge consisitng of armchair and zigzag segments 
whose typical length dg significantly exceeds the elec- 
tronic wavelength, de S> Xe, where Ac = 2v/ujin- This 
is an intermediate case between the two limiting cases 
considered in the previous two sections. 

Only armchair segments contribute to the Raman 
process. Moreover, contributions from different seg- 
ments add up incoherently. Thus, we first focus on 
the contribution of a single armchair segment, placed at 
a; = 0, —de/2 < y < de/2 (as before, graphene is as- 
sumed to occupy the half-space a; > 0). The electronic 
Green's function corresponding to the reflection from a 
single armchair segment can be easily written from the 
Huygens-Fresnel principle, Eq. (j3.7p if one approximates 
its value on the boundary by that for an infinite perfect 
armchair edge: 



de/2 



G(r,r';e) = Go(r - r',e) - J dye 



-de/2 



X Go{x, y - ye; e) Y.^ Go{x', ye - y'] e) S^Ay. 

(8.1) 



This approximation ignores the change of the exact wave 
function within the distance ~ Ac from the ends of the 
segment, which gives a small error if c?e ^ Ac. In fact, it 
is the standard approximation for the study of diffraction 



m optics 
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Using this Green's function, we obtain the following 
expression for the matrix element corresponding to emis- 
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sion of a phonon in an arbitrary point (X, Y) : 



M = - \ 




dye 



Pin P out P[„ Pout 



X COS ■ 



i^out 

sm 



Pa 
P'a 



^{X - XaY + {ya - YY 

ya-ye , 



arctan ■ 



arctan 



yg-Y 

X ~Xa'' 
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■Pa 



■ 27 



(8.2a) 

(8.2b) 
(8.2c) 

(8. 2d) 
(8.2e) 



V V 

where a is either "zn" or "ouf . It is convenient to use the 
paraxial approximation with respect to the axis connect- 
ing the points (0, ?;e) and {X,Y). In this approximation 
we expand 

XcOS^(j)0 iVa-yaoY 



Pa+ Pa 



X 



cos 00 Xa{X - Xa) 

ya ~yaO 2, 

cos 00, 



Va - VaO 



" X-Xa 

yao = tan0o, tan^o 



cos 00, 

Y-ye 



X 



. -Sa) 
(8.3b) 
(8.3c) 
(8.3d) 



Integrating over Tin , ro„t in the usual way, we obtain 



cjphX/(2w COS 0o) 



X exp 



^out 



27 



X 



2 / u cos 00 

X s\n{ipin - 0o) sin((^out - 0o)- 



(8.4) 



This integral can be calculated analogously to the 
standard diffraction integral in optics4i According to 
Eq. ()5.10a|) . the integrated intensity into the full solid 
angle 47r, summed over two polarizations of the emitted 
photon, and integrated over R, is given by 



Id 



8\k /e2\%2 



2 T T 2 



X < siiY ipin In ■ 



Ul. 



ph 



(47) 



(47)^ ^ , 

■ + COS (Pi; 



^In^ 



2v 



■In 



(47)^ 



^^ph 



(47)2 



(8.5a) 



This expression is analogous to Eq. (j6.5bp . weighted 
by the factor de/Ly. The coefhcient at the last term, 
cx cos^ ifiin , determines the minimum of the intensity at 
ipin = and has two contributions: the one correspond- 
ing to the infinite edge, and the one due to the finite 
size of the segment. The latter one is dominant unless 
de ^ (w/win)('^m/wph)2 ~ 50 nm for tOin = 2eV. Still, as 
long as LOinde/v ^ 1, the ratio between the intensities for 
the parallel and perpendicular polarizations is large. If 
we use Eq. (|5.10b|) for the intensity in a solid angle Oout 
in the presence of an analyzer, we obtain 



= 4A 



9 \ 2 9 

Oout I e \ w 
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An 



LUinLx Ly Wpjj 



X ■! sin ipin sin ipout In 



■^ph 



(47)2 



(47)2 



+ 



siT?{(pin + ipout) + 2 sil'i sin 2ipout 
V , ^ph + (47)' , 



■In- 



(47)2 



-t- COs2 C0S2 Lpout 



ph j^^^ ^in 



" ■In'P^ 



(47)2 



(47)2 



(8.5b) 



Eqs. (|8.5ap . (|8.5b[) describe the contribution of a single 
armchair segment to the Raman intensity. To obtain the 
contribution of the whole edge, it is sufhcient to multiply 
these expression by the total number of such segments 
and replace dg by its average, if all segments have the 
same orientation. It is crucial, however, that up to three 
different orientations of armchair segments are possible, 
at the angle 7r/3 to each other, as discussed by the author 
and coworkers in Ref . [l9l 

Let us first consider a measurement in the absence of 
an analyzer. Since the intensity is a biliniear form of the 
polarization vector e^n, it can always be written in the 
form 

I{ipin) OC C0s2((^i„ - ip-max) +eSlV?{Lpin ~ Ifimax) 

= + ~ 2l^™aa;), (8.6) 

where (pmax is the angle where the intensity is maximum 
and £ is the ratio between the intensities in the minimum 
and in the maximum. Eq. (j8.5ap corresponds to fmax ~ 
7r/2 and small e ^ 1 due to the quantum diffraction. 

Let the edge have iVo armchair segments oriented along 
the y direction, such as the one considered above, and N± 
segments oriented at ±7r/3 to the y axis. Note that the 
average direction of the edge may still be arbitrary, as 
it depends on the distribution of zigzag segments too. 
Let each segment be charaterized by the same values of 
e and (pmax^ when the latter is measured with respect to 
the corresponding normal. Adding the contributions, we 
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again obtain an expression of the form of Eq. (|8.6p , but 
with different values of parameters: 

NoI{ip^n) + N+I{ip,n ~ 7r/3) + iV_/((p,„ + 7r/3) 

OC COS^(<^i„ - ipmax) + SSm^{ipin ~ V>max), (8.7a) 

-arctan ^ + '-, 8.7b) 

2 2No-N+-N^' ^ ' 



(l+£)jVtot-(l-£)iV 
(l + e)Artot + (l-£)iV' 



(8.7c) 



N = ^JNl^ - i{N+N^ + iVoiV+ + 7VoiV_), (8.7d) 
A^toi = A^o + A^+ + A^-. (8.7e) 



Inspection of Eq. l|8.7bp shows that e ^ 1 if and only if 
(i) e < 1 and (ii) Ntot - < Ntot- The latter condition 
is equivalent to having one of Nq, N+, N- to be much 
larger than the others. If these conditions hold, we can 
write (assuming that iVo ^ , iV_ for concreteness) 



3 iV+ + iV_ 



Nn 



In the opposite case Nq 



iV„ we have e = I, 



so that isotropy is fully restored and no signatures of 
quantum diffraction are left. 

An analogous summation can be performed in the pres- 
ence of an analyzer in the general case, but the final ex- 
pressions are very bulky and not very informative. The 
qualitative conclusion is the same: the terms which were 
small compared to the leading term sin'^ (fin sin^ ipout 
grow as one adds segments with different orientations. 
At A^o = A^+ = A^- the isotropy is restored, 

sin^ ipin sin^ ipout + sin^(v3i„ - 7r/3) sm^{ipout - 'I'/S) 
+ sin^((^i„ + 7r/3) sm^{ipout + 7i'/3) 

3 9 

= - sin^(i^i„ - ipout) + ^ cos^((pi„ - (font), (8.9) 
o o 

and signatures of the quantum diffraction are lost. 

Let us focus on the special case when the average direc- 
tion of edge is zigzag, and it is symmetric on the average, 
A^+ = A^o, ^- = 0. Then we have 

Id cx sin^ ipin sin^ ipout 

+ sin^ (^(pin - ^) sin^ (^ipout ^ ^) + ^i^) 
3 3 

^8^4 ~ 

cos ( - g ) cos 



(8.10) 



The maximum of intensity is reached when both polar- 
izations are along the average direction of edge. For 
the unpolarized detection, we add the contributions with 
font = and (pout = ''■/2 and obtain the ratio be- 
tween the minimum and the maximum intensity to be 
1/3 -I- 0{e); if an analyzer is used and ipin = ipout, we 



obtain e = 1/9-1- 0(e). These findings agree with the 
available experimental daXaM^ 

The dependence of Eq. (|8.10p has a remarkable prop- 
erty: at (fiin = 0, (fiout = T^l'i (or vice versa) the leading 
term vanishes and lu = 0(e), i. e. the quantum limit 
is still accessible. In fact, the same will be true for any 
edge with only two orientations of the segments (i. e, for 
A^_ — 0, but A^o ^ A^+, generally speaking). The ratio 
of intensity in this minimum to the maximum without 
analyzer is given by 



ph_ ^ph 



In 



^ph 



(47)' 



(47)2 



^Nl + Ar2 - NqN^ 



.11a) 
.lib) 



IX. CONCLUSIONS 

We have studied scattering of Dirac electrons on a 
graphene edge. For a translationally invariant edge (such 
as zigzag or armchair or another edge with a certain spa- 
tial period) the reflection can be described by an effective 
low-energy boundary condition for the electronic wave 
function. 2^-2^ For edges which are rough on the atomic 
scale we have proposed a random-matrix model which 
describes random scattering of electrons on the edge, re- 
specting the particle conservation and time-reversal sym- 
metry. Essentially, each point of the edge acts as an in- 
dependent point scatterer with a random rotation of the 
valley state. We have also considered edges consisting 
of zigzag and armchair distinct segments longer than the 
electron wavelength, each segment can be treated as a 
piece of an ideal edge, while the small corrections due 
to quantum diffraction, can be found using the Huygens- 
Fresnel principle for Dirac electrons, analogously to the 
standard treatment of diffraction in the classical optics. 

Next, we have calculated the intensity of the edge- 
induced D peak in the Raman scattering spectrum of 
graphene. It is shown how the quasiclassical character 
of the electron motion manifests itself in the polarization 
dependence of the intensity. For an ideal armchair edge 
the maximum of intensity corresponds to the case when 
both the polarizer and the analyzer are along the edge, 
and the large ratio of intensities in the maximum and 
the minimum turns out to be determined by the quan- 
tum corrections to the quasiclassical motion of the pho- 
toexcited electron and the hole. For an edge consisting 
of randomly distributed zigzag and armchair segments 
of the length significantly exceeding the electron wave- 
length, the effect of quantum diffraction can be masked 
by the presence of armchair segments of different orienta- 
tions. The maximum and the minimum of the intensity 
are determined by the number of the armchair segments 
with different orientations, rather than the average di- 
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rection of the edge. If only two orientations of armchair 
segments are present in the edge, the quantum diffraction 
hmit can still be probed by a careful choce of the polar- 
izer and the analyzer (the polarizer should be oriented 
perpendicularly to one of the armchair directions, the 
analyzer perpendicularly to the other one). For an edge, 
rough at the atomic scale, no segments can be identified, 
and the intensity reaches its maximum for the polariza- 
tion along the average direction of the edge. The ratio of 
the maximum and the minimum intensity is determined 
by the ability of the edge to scatter electrons at oblique 
angles. 

As the whole Raman process is edge-assisted, one can 
pose the question about the characteristic length scale 
which restricts the process to the vicinity of the edge. 
We find that the answer is not unique, and the length 
scale depends on the specific observable under study. If 
one is interested in the total intensity or its polarization 
dependence, the effective length scale is w/wph {v being 
the electronic velocity and cjph the phonon frequency). 
However, if one makes a spatially resolved experiment, 
measuring the dependence of the intensity on the position 
of the excitation spot, the relevant length scale is the 
electron inelastic scattering length v/{2^). We have thus 
found a qualitative agreement with the interpretation of 
Ref. IG, but we argued the inelastic scattering length 
found in that work is too large to be consistent with other 
available information on electron inelastic scattering in 
graphene. 



X. ACKNOWLEDGEMENTS 



correlated with any cristallographic direction, so the de- 
scription is equivalent for all orientations. Thus, expres- 
sions for a general orientation of the edge are obtained 
by replacing the coordinates by x ^ n • r, y ^ [n x r]^, 
and the polar angles by if —^ ip + ip^. 

Let us label the eigenstates of the problem by three 
quantum numbers: (i) e, the energy of the electron; 
(ii) Py, the y-component of the momentum of the inci- 
dent plane wave (obviously, \py\ < |e|/f), (iii) k, the val- 
ley index (we write «; = ±1 for AT and AT', respectively). 
The most general wave function of such an eigenstate has 
the form 



V'£,p„,K(a;,y) 





V Wx 

p,^^{e/vr-p2, p'^ ^ ^{e/vY - [p'yY, (A2b) 
_ f Lydpy „ 27r 
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Sipy-p') (A2c) 



-\e\/v 



sgn e 



'4' <i,Py,K — Sy'0e,Py,K 



1 
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e ^'^pa/^sgne 



arctan 



Py 



2 - '^Py - 2' 



(A2d) 

S)((^e) 
(A2f) 
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APPENDIX A: SCATTERING MATRIX FOR AN 
EDGE 

The representation of the 4-column vector ip, natural 
for scattering problems, is that of Ref. IH, Eq. (|3.3p . For 
4-column vectors which can be represented as a direct 
product 



xiyi 



X2y2 



X2yi 








yi 




Xl 


xiy2 




X2 




y2 




X2 



® iyiff'K + y2'f>K' 



(Al) 

the S matrices act on the x variables, while the A ma- 
trices act on the y variables. The basis in the valley 
subspace is denoted by 4>k, 4>k' for future convenience. 

To keep the formulas compact we assume that the edge 
is along the y direction and graphene is occupying the 
half-space a; > 0, i. e., n = e^;. The average orientation 
of a disordered edge with de ^ a does not have to be 



The prefactor {vpx) in Eq. (|A2ap fixes the com- 

ponent of the probability current perpendicular to the 
boundary to be jx = —vsgTie/{LxLy) for the incident 
wavei^ Finally, S3 !^ are the scattering matrix elements, 

PyPy 

which may depend on e. Note that ipp — ippi corre- 
sponds to the specular reflection. The wave functions 
i'e,py,Kix,y) are orthogonal. 



oo oo 



dy I dxiplp^^^{x,y)ilJe',p'y,K'{x,y) = 

5{e-e')dp^p,5^,,, (A3) 



-oo 
2'KV 

Lx 



provided that the scattering matrix is unitary: 

E {s;(p^*s;(^: = s..„6p^p.. (A4) 

P'y,K' 

The same condition ensures the conservation of the flux, 
i. e. that the current = vsgne/ (L^Ly) for the scat- 
tered part of the wave function (jA2ap . Reflection from 
a regular edge with wave function given by Eq. 

corresponds to S"^/^^ = <5p^Py'/'I"S'A0K. 
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The time reversal symmetry imposes another condition 
on the scattering matrix (reciprocity condition) . Namely, 
the time-reversed wave function Utipl .p^^^{x,y) describes 
an eigenstate of the problem, and thus it must be a linear 
combination of ip^p' ,k' (a;, y) with different p'y, n' . Noting 
that 

C/tV'e -Py -K = ~-ii^i^e,py,K- (A5b) 
we obtain the reciprocity condition. 



last term of Eq. (|4.6p , we obtain the relation 



P'yPv 



--SZ^p",,, (A6) 



which can also be written in the matrix form as Sn' n = 
—AyS'^„ „, Ay, where S*^ denotes the 2x2 matrix 

^ Py^Py ^ 

transpose. 

To relate the scattering matrix to the T-matrix in 
Eq. (|4.12p . we compare the x ^ ti/|e| asymptotics of 
the corresponding wave functions. The Green's function 
Go(r — re,e) entering Eq. I|4.12p can be represented as 
follows: 



iPy{y-Vc) 



Go{x,y - ye,e) = 

oo W\/v 

de' r dpy\£\_ 

27rv J 27r vpx e — e' + zO+ sgn e 

-\e'\/v 

K 

(A7) 

where we suppressed the damping 7^ — > 0+ and denoted 
Px = \J (e'/u)2 — p^) . To determine the x ^ v/\t \ asymp- 
totics of the Green's function, we shift the e' integration 
contour in the upper (lower) complex half-plane at e' > 
(e' < 0) for the term cx 6*^="^, and in the lower (upper) 
half-plane at e' > (e' < 0) for the term cx e^*P-"=. The 
magnitude e" of the shift is such that v/x <C e" ^ |e'| 
As a result, (i) the contribution from the term oc e~*P"=^ 
will be exponentially small except for the small region 
near zero, |e'| ~ v/x; (ii) the contribution from the term 
cx e*P^^ will be determined by the pole e' = e + iO"*" sgn e, 
the rest of the contour contributing an exponentially 
small quantity; (iii) the contribution from |e'| ~ v/x will 
be small as v/{\e\x). Thus, at x » w/|e| we arrive at 



\e\/v 



Go{x,y- ye,e) 



f,iPyiv-y'!)+W=.x 



V J 2lT VPx 



(A8) 



Substituting this expression and the first term of 
Eq. (|4.6p into Eq. (|4.15p . and comparing the result with 



K K 

P'yPy 



^^g"^ f ^ ^-^^p'y-Py)y. 



'^l,p'y,K'T{^p'y , - (^p„ ; ye)^e,py,^ 



.(A9) 



Here ippi is the polar angle of the direction s, and ir — ipp^ 
is the polar angle of the incident momentum p, directed 
towards the edge. Using Eqs. (|A5ap . (|A5bp . one can es- 
tablish the equivalence between the reciprocity condition, 
Eq. (jA6p . and Eq. (|4.13p . Using the unitarity condition, 
Eq. (jA4p . one fixes the prefactor Trw/jel in front of the 
(5- function in Eq. (|4.14dp . 



APPENDIX B: INTEGRALS OF SEC. [yTBl 



First, let us focus on Eq. (|6.3ap for Upon inte- 

gration over yi„, yout, it takes the following form: 



Mn = — 



00 X 



4Ar<-u f de 



-00 

•^in ~^ ^out ^ A^ ^in ^ A Xout 

|e| \(^in - e| \^^out - e| 

9x ^ — {Xin + Xout ) H (A - X,,n ) 

V V 



(A - Xout)- 



-1/2 



(Bib) 



If one neglects the square root and integrates first over 
Xin, the energy is constrained to be e « with 
the precision ~ v/X. If one integrates over Xout, then 
e w LUout/'2 (which is not inconsistent with the previous 
condition as long as A ~ w/cjph). If one also removes the 
moduli and integrates over e, then Xin + Xout ~ A. 

To implement these observations more rigorously, let 
us introduce new variables, 

~ -^in ~^ -^out A /Tiio ^ 

X = , Xf) = Xin Xout J [tjZa.) 



(B2b) 



V 2v 



2v 



Recalling that the typical scale of the e dependence of the 
rest of the integrand is u>in, and that of the x dependence 
is X, we rewrite 



ex LOinX ( e 



X \ LOinX I e X 



V \ LOin X 



V \ UJin X 



m 

Taking the two expressions in the brackets as new inte- 
gration variables, we can use the stationary phase ap- 
proximation justified by tOinX/v 3> 1. As a result, the 
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integral is contributed by |e| \/ uJinV j X <^ ujin and 
by \x\ ^ \J XvjijJin <C X. Thus, in the rest of the inte- 
grand they can be simply set to zero, i. e. we approximate 
g4je2;/« ^ (7rw/2)(5(e)(5(i). The remaining integration over 
xq is elementary: 



X {X~-\xa\)l2 

dxo I dx 



di 



-X -(X-\xo\)/2 \e\<^,,j2 



For the integral in Eqs. (|6.6ap . (j6.6b[) we have: 




FIG. 8: (Color on-line.) Geometry corresponding to Eq. (|B7 
The solid black arrows show the vectors pi,p2,P3. 



~ •^in ^out •^in ~t~ •^out ^/T-»r- ^ 

X = , xo = ^ (,B5a) 



€ — € 



Win -|- UJout 



^ -r -r 

V V V 



(B5b) 
(B5c) 



problem] , 



M 



— d r^nd Tout COS ■ 

V . 2 



dl 



X i/2 



— J dx J dxQ e 

|e|<aj„/2 -x/2 



^-i(2e-i^ph/2)X/v _ I 

21 - a;ph/2 



— (^ph — ^ph) 



(B6) 



If we extend the limits of the e integration to infinity, the 
integral vanishes by analyticity. Thus, it is determined 
by the region |e| ^ Wi„ and can be estimated as ~ v/ujin- 

Now let us pass to the derivation of Eq. (j6.4bp for 
M±. The difference from the previous case is that in- 
tegration over ijin , Hout in the derivation is more sub- 
tle. Namely, we encounter oscillating integrals of the 
kind / j/^e*"** dy. This integral is assigned the value 
^/^^fA{—ia)~^/'^, which may be understood as the an- 
alytical continuation from the upper complex half-plane 
of a. One may argue that since the integral is divergent 
for real a, it is determined by \y\ ~ 1 / \/Im a which are 
not small, so the expansion to both in the exponential 
and in the pre-exponential factor is not valid. 

Let us write down the expression for the matrix ele- 
ment without using the paraxial approximation [we omit 
the 1 /C terms of the expansion in Eq. p.6p as they do not 
produce y'^ terms and do not represent any convergence 



. 2ipin + >^l- ^2 . 2ipout - ^1 - ^3 

X sm sm 

2 2 

V — l)pi/'V+{i\uJi„-<i\ — i)p2/v + {i\LJout—e\—'f)p3/v 



{ujin - e){uJout - e.)e\ 

P1P2P3 



(B7) 



The distances jOi,2,3 and the angles (pi.2,3 are shown 
in Fig. [51 Again, integration over e/v gives effectively 
5(pi — P2 — Pz), so the exponential in the integrand be- 
comes e(*'^".~27)P2/i'-H(«u;„„t-27)P3/»', The phase of this ex- 
ponential is rapidly oscillating, and the stationary phase 
condition gives precisely yimVout — > 0. The straight- 
forward application of the stationary phase method for 
^in = font = is impeded by the presence of two sine 
functions in front of the exponential. 

First of all, let us check the convergence of the integral 
at large distances, i. e. assuming Tin, rout S> X. Then 
it is convenient to use the polar coordinates for ri„ , Tout 
whose polar angles are almost tt + (p2,TT + tpz- The S- 
function fixes ip2 = tt — ipi, (^3 = tt + </ji . We can set 
X = {) everywhere except the Jacobian: 



Pi 



1 '^out 



'^in ^out 



P2,3 ~ rin^out + XC0S(P2,3, 



SiPl - P2- P3) 



{ip2 + (^3)^ (B8a) 
(B8b) 
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This gives a finite value, 

7r/2 



M 



-7r/2 



d<^i / nndrin rout dr out 



(^,„/2)3 



y g(iWi„-27)ri„/t)+(j:i^o„t-27)rout/t) 

X C0S(/5i COs(<^i„ + COs((/Jo„4 - (fi) 



87r3 Wj„X Win 



X suiipinSaiipout - ^ cos(^i„cosv9o„t^B9) 

Note that for any polarization it is smaller than Ai\\ from 
Eq. (|6.4ap by a factor Wph/wi„. 

Thus, we are facing a situation of the kind 



(BIO) 



where (i) the function $(y) is growing as |y| at |y| ^ oo 
and has an extremum at y = 0, where it can be ex- 
panded as ^{y) = $(0) + $"(0) 2/2/2 + 0{y^), (ii) a is 
large and is assumed to have a positive imaginary part, 
and (iii) f{y) = f"{0)y^/2 + 0[y^) at y ^ 0, while the 
behavior of the function /(y) at |y| ^ oo is such that the 
in the limit Ima — > the integral F(a) remains finite. 
Then F{a) is analytic in the upper complex half-plane 
of a. Thus, the expansion 



F{a) 



A72r(o) 

-za$"(0)]3/2 



o{{- 



^-5/2^ 



(Bll) 



established for large positive imaginary a, holds in the 
whole upper half-plane, including the vicinity of the real 
axis. The two-dimensional integration over yin,yout, 
leading to Eq. (|6.4bp is fully analogous. 



APPENDIX C: DETAIL ED EX PRESSIONS FOR 
SEC. IVICI 

Note that length I2 is not well defined in the general 
case. Namely, for an arbitrary kernel IC{x,x') and an ar- 
bitrary shape of Sin {x) the expansion in the derivatives 
of £in{x) does not automatically "wrap" into the expan- 
sion in the derivatives of ^(a;), Eq. (|6.9p . However, I2 can 



be defined in two particular cases which are most impor- 
tant for us. It is convenient to introduce an auxiliary 
notation: for an arbitrary function T{x, x') we define 



/fix rlx^ 
T(x,x')K{x,x') - HhM)- 



(CI) 

One case when the length I2 can be defined is when the 
kernel IC{x,x') is such that 



{{x')) + {{{x'r)) = 'iiixx')), 



(C2) 



then one can set I2 = {{xx')). For the kernel given by 
Eq. (|6.7cp this property holds in the limit £^ ^ ^ph: 



x'^]C{x, x') dx dx' 



e-2 



, p2p2 7 



ph + 20^2,2 



7 ph 



4(^^+4£2J2 

00 

j xx' lC{x ^ x' ) dx dx' = if^Yi^Ji 



(C3) 



p2 p2 f2 
^7 



02 
•^ph 



(C4) 



The second case when I2 can be defined is the gaussian 
profile fin (a;) cx e^(^^^°) /'^'^^ \ which has a special prop- 
erty £'^^{x)£M = KJa;)]^ - [£,n{x)f / L\ Then the 
difference between the left-hand and the right-hand sides 
of Eq. (|C2p can be absorbed in the overall coefficient, so 
that instead of Eq. (|6.7cp we have its slightly modified 
version: 



w(li 



w"{li 
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{{x^)) + {{{x'f))~2{{xx')) 
4L2 



(C5a) 
(C5b) 
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{{x^)) + {{{xr))+2{{xx')) 



^2(^4 + 10^2^2^ + 8^4,)/ £2,^ 4^: 



4(^^ + 4^^h)' 
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^1 



4(^2 4^2 J2 
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For the states with e < the direction of the electron 
group velocity is opposite to that of momentum. Thus, for 
e < the first term of the wave function (|A2a|) represents 
an outgoing flux of electrons, and the second term - to 
an incoming flux. However, we prefer to work with the 
chronologically ordered Green's function G(p,e), defined 
in Eq. (|3.4p . rather than for the retarded Green's function 
G'^(p,£), the analytical continuation of G'(p, e) from the 
positive semiaxis of e through the upper complex half-plane 
of e. For this purpose we work with the solutions (IA2a|) 
where the first term of the wave function (IA2a|) represents 
an incoming fiux for electrons if e > and for holes if e < 0. 



